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ABSTRACT 

This  paper  reports  the  outcome  of  an  exhaustive  analyt¬ 
ical  and  numerical  Investigation  of  stability  and  ro¬ 
bustness  properties  of  a  wide  class  of  adaptive  control 
algorithms  in  the  presence  of  unmodeled  dynamics  and 
outpuc  disturbances.  The  class  of  adaptive  algorithms 
considered  are  those  commonly  referred  Co  as  model- 
reference  adaptive  control  algorithms,  self-tuning 
controllers,  and  dead-beat  adaptive  controllers;  they 
have  been  developed  for  both  continuous- time  systems 
and  discrete-time  systems.  The  existing  adaptive  con¬ 
trol  algorithms  have  been  proven  to  be  globally  assymp- 
tocically  stable  under  certain  assumptions,  Che  key  ones 
being  (a)  that  the  number  of  poles  and  zeroes  of  the 
unknown  plane  are  known,  and  (b)  that  the  primary  per¬ 
formance  criterion  is  related  to  good  command  following. 
These  theoretical  assumptions  are  coo  restrictive  from 
an  engineering  point  of  view.  Real  plants  always  con¬ 
tain  unmodeled  high- frequency  dynamics  and  small  delays, 
and  hence  no  upper  bound  on  the  number  of  the  plant 
poles  and  zeroes  exists.  Also  real  plants  are  always 
subject  to  unmeasurab le  output  additive  disturbances, 
although  these  may  be  quite  small.  Hence,  it  is  impor¬ 
tant  to  critically  examine  the  stability  robustness 
properties  of  the  existing  adaptive  algorithms  when 
some  of  Che  theoretical  assumptions  are  removed;  in 
particular,  their  stability  and  performance  properties 
in  the  presence  of  unmodeled  dynamics  and  output  dis¬ 
turbances  . 

A  unified  analytical  approach  has  been  developed  and 
documented  in  the  recently  completed  Ph.D.  thesis  by 
Rohrs  CIS ]  that  can  be  used  to  examine  the  class  of 
existing  adaptive  algorithms.  It  was  discovered  that 
all  existing  algorithms  contain  an  infinite-gain  opera¬ 
tor  in  the  dynamic  system  that  defines  command  ref¬ 
erence  errors  and  parameter  errors;  it  is  argued  that 
such  an  infinite  gain  operator  appears  to  be  generic  ts 
all  adaptive  algorithms,  whether  they  exhibit  explicit 
or  implicit  parameter  identification.  The  phacXLtat 
engineering  con sequences  of>  the  existence  o&  the  in- 
AAotiXe-gain  opvuUox  axe.  diiOAhoui.  Analytical  and 
simulation  results  demonstrate  that  sinusoidal  referenoe 
inputs  at  specific  frequencies  and/or  sinusoidal  output 
disturbances  at  any  frequency  (including  d.c.)  cause  the 
loop  gain  of  the  adaptive  control  system  to  increase 
without  bound,  thereby  exciting  the  (unaodeled)  plant 
dynamics,  and  yielding  an  unstable  control  system. 

Hence,  it  is  concluded  that  none  of  the  adaptive  algo¬ 
rithms  considered  can  be  used  with  confidence  in  a 
pnacXical  control  system  design,  because  initability 
utill  set  in  i sith  a  high  pnobability. 
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1.  INTRODUCTION 

Due  to  space  limitations  we  cannot  possibly  provide  in 
this  paper  analytical  and  simulation  evidence  of  all 
conclusions  outlined  in  the  abstract.  Rather,  we 
summarize  the  basic  approach  only  for  a  single  class 
of  continuous-time  algorithms  that  Include  those  of 
Monopoli  [4],  Narendra  and  Valavanl  [1],  and  Feuer  and 
Morse  [2].  However,  the  same  analysis  techniques  have 
been  used  to  analyze  more  complex  classes  of  (1)  con- 
tlnuous-tima  adaptive  control  algorithms  due  to 
Narendra,  Lin,  and  Valavanl  [3],  both  algorithms  sug¬ 
gested  by  Morse  [4],  and  the  algorithms  suggested  by 
Egardt  [7]  which  include  those  of  Landau  and  Silveira 
[6],  and  Kreisselmeler  [19];  and  (2)  discrete-time 
'—adaptive  control  algorithms  due  to  Narendra  and  Lin  [22], 
Goodwin,  Ramadge,  and  Caines  [23]  (the  so-called  dead¬ 
beat  controllers),  and  those  developed  in  Egardt  [17], 
which  Include  the  self-tunning-regulator  of  Astrom  and 
Wittenmark  [18]  and  that  due  to  Landau  [20],  The 
thesis  by  Rohrs  [15]  contains  the  full  analysis  and 
simulation  results  for  the  above  classes  of  existing 
adaptive  algorithms. ' 

The  end  of  the  1970’s  marked  significant  progress  in 
Che  theory  of  adaptive  control,  both  in  tens s  of  ob¬ 
taining  global  asymptotic  stability  proofs  [1-7]  as 
well  as  in  unifying  diverse  adapeive  algorithms  Che 
derivation  of  which  was  based  on  different  philosophical 
viewpoints  [8,9]. 

Unfortunately,  the  stability  proofs  of  all  these  algo¬ 
rithms  have  in  common  a  very  restrictive  assumption. 

For  continuous-time  Implementation’ this  assumption  is 
that  Che  number  of  the  poles  and  zeroes  of  che  plane, 
and  hence  its  relative  degree,  i.e.,  its  number  of  poles 
minus  its  number  of  zeroes,  is  known.  The  counterpart 
of  this  assumption  for  discrete-time  systems  is  that 
che  pure  delay  in  the  plant  is  exactly  an  integer  number 
of  sampling  periods  and  that  this  integer  is  known. 

This  restrictive  assumption,  in  turn,  is  equivalent  to 
enabling  che  designer  to  realize  for  an  adaptive 
algorithm,  a  positive  real  error  transfer  function,  on 
vhich  all  stability  proofs  have  heavily  hinged  to-date 
[8J.  Positive  realness  Implies  Chat  che  phase  of  the 
system  cannoc  exceed  +  90*  for  all  frequencies,  while 
it  is  a  well  known  fact  chat  models  of  physical  systems 
become  very  Inaccurate  in  describing  actual  plant  high- 
frequency  phase  characteristics.  Moreover,  for  prac¬ 
tical  reasons,  most  controller  designs  need  to  be  based 
on  models  which  do  not  contain  all  of  che  plane  dynamics, 
in  order  to  keep  the  complexity  of  che  required  adaptive 
compensator  within  bounds. 

Motivated  from  such  considerations,  researchers  in  che 
field  recently  started  investigating  the  robustness  of 
adaptive  algorithms  to  violation  of  the  restrictive 
(and  unrealistic) assumption  of  exact  knowledge  of  the 
plant  order  and  its  relative  degree.  loannou  and 
Kokotovlc  [10]  obtained  error  bounds  for  adaptive 
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observers  and  identifiers  la  che  praaanca  of  unmodeled 
dynamics,  whlla  auch  analytical,  results  wara  hardar  co 
obtain  for  raducad  ordar  adapt lva  coatrollari.  Tha 
first  auch  result,  obtalnad  by  Rohrs  at  al  [11] , 
consist a  of  "linearization"  of  tha  arror  equations, 
undar  cha  assumption  that  tha  ovarall  system  la  In  lta 
final  approach  to  convergence.  loannou  and  Kokotovic 
[ 12]  latar  obtalnad  local  stability  raaulta  In  tha 
praaanca  of  unmodalad  dynamics ,  and  should  that  cha 
apaad  ratio  of  alow  versus  fast  (unaodalad)  dynamics 
directly  affected  cha  stability  region.  Earliar  simu¬ 
lation  studios  by  Rohrs  at  al  (13}  had  already  shown 
increased  sensitivity  of  adaptive  algorithms  to  dis¬ 
turbances  and  unaodalad  dynamics,  generation  of  high 
frequency  control  Inputs  and  ultimately  Instability. 
Simple  root- locus  type  plots  for  tha  linearized  system 
In  [11]  showed  how  tha  presence  of  unaodalad  dynamics 
could  bring  about  Instability  of  tha  ovarall  system-  it 
was  also  shown  there  chat  tha  generatsd  frequencies  In 
cha  adaptive  loop  depended  nooliaaarly  on  tha  magnitudes 
of  tha  reference  lnpuc  and  output. 

Tha  mein  contribution  of  this  paper  is  in  showing  that 
two  operators  inherently  included  in  all  algorithms 
considered  —  as  part  of  the  adaptation  mechanism  — 
base  Infinite  gain.  As  a  result,  two  possible  mech¬ 
anisms  of  instability  are  isolated  and  discussed.  It  la 
argued,  chat  the  destabilizing  effects  in  the  presence 
of  unaodalad  dynaadcs  can  be  attributed  to  either  phase 
—  in  the  case  of  high  frequency  inputs  —  or  primarily 
gain  considerations  —  In  the  case  of  unmeasurable  out¬ 
put  disturbances  of  any  frequency.  Including  d.c.,  which 
result  in  nonzero  steady-state  errors.  The  latter  fact 
is  most  disconcerting  for  Che  performance  of  adaptive 
algorithms  since  it  cannot  be  dealt  with,  given  that  a 
persistant  disturbance  of  arty  frequency  can  have  a  des¬ 
tabilizing  affect. 

Our  conclusions  that  che  adaptive  algorithms  considered 
cannot  ba  used  for  practical  adaptive  control,  because 
t he  physical  system  hri.ll  eventually  become  unstable , 
are  baaad  upon  two  facts  of  life  that  cannot  be  ignored 
in  any  physical  control  design:  (1)  there  are  always 
unaodeled  dynamics  at  sufficiently  high  frequencies 
(and  it  is  futile  to  cry  to  model  unaodeled  dynamics) 
and  (2)  the  plant  cannot  be  isolated  from  unknown  dis¬ 
turbances  (a.g.,  60  Hz  hum)  even  though  these  may  be 
small.  Neither  of  these  two  practical  isauas  have  been 
Included  in  tha  theoretical  aasumptlons  common  to  all 
adaptive  algorithms  considered,  and  chis  is  why  these 
algorithms  cannot  be  used  with  confidence.  To  avoid 
exciting  unaodeled  dynamics,  stringent  requirements 
must  be  placed  upon  the  bandwidth  and  phase  margin  of 
cha  control  loop;  no  such  considerations  have  been 
discussed  in  the  literature.  It  is  not  at  all  obvious, 
nor  easy,  how  co  modify  or  extend  the  available  algo¬ 
rithms  to  control  their  bandwidth,  auch  less  chair 
phase  margin  properties. 

In  Section  2  of  chis  paper  proofs  for  Che  Infinite  gain 
of  the  operators  generic  to  the  adaptation  mechanism 
are  given.  Section  3  contains  che  development  of  two 
poaelble  mechanisms  for  instability  chat  arise  as  a 
result  of  tha  infinite  gain  operators.  Simulation 
results  chat  show  the  validity  of  che  heuristic  argu¬ 
ments  in  Section  3  are  presented  in  Section  A.  Section 
5  contains  cha  conclusions . 

2.  THE  lUDl  MODEL  STRUCTURE  FOR  A  REPRESENTATIVE 
ADAPTIVE  ALGORITHM 

The  simplest  prototype  for  a  model  reference  adaptive 
control  algorithm  in  continuous- tine  has  its  origins  to 
at  least  as  far  back  as  1974,  la  tha  paper  by  Monopoli 
(14],  This  algorithm  has  been  proven  asymptotically 
stable  only  for  the  case  when  the  relative  degree  of 


che  plant  Is  unity  or  at  most  two.  The  algorithms 
published  by  Karendra  and  Velavani  [1]  aid  Fetter  and 
Morse  [2]  reduce  to  the  same  algorithm  for  che  perti¬ 
nent  case.  This  algorithm  will  henceforth  be  referred 
to  as  CA1  (continuous-time  algorithm  No.l) 


The  following  equations  summarize  the  dynamical  equa¬ 
tions  that  describe  it;  see  also  Figure  1.  Tha  equa¬ 
tions  presented  here  pertain  co  the  case  where  a  unity 
relative  degree  has  been  normally  assumed.  In  the  equa¬ 
tions  below  r(t)  is  the  (command)  reference  input,  end 
d(t)-0. 

Plant: 


y(t) 


gpB(a) 

A(s) 


(u(t) ] 


(1) 


Variables :  «ui(t>-pl7r  . -l  (2) 
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Equation:  e(t)  • 
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In  Che  above  equations  the  following  definitions  apply: 


k(t)  =  k*  +  E(t) 


(10) 


where  k*  is  a  constant  2n  vector 


K*(s) 


kS(n-l/'2  +  k;(n-2/_3+-*Sl 


where  k*^  is  the  1-ch  component  of  k* 


K*(s) »  k*  a"*1  + 

y  y« 


’'y(n-l) 


where  k*^  is  the  1-th  component  of  k*  and  the  vector 

k*  componenwise  corresponds  exactly  to  the  vector  k(t) 
In  eqn.  (3a).  In  the  preceding  equation  we  have  tried 
to  preserve  the  conventional  literature  notation  (3,4, 
5,9],  with  P  representing  the  characteristic  polynomial 
for  the  stats  variable  filters  and  k(t)  the  parameter 


misalignment  vector. 


The  quantity  a*B*  rapraaants  the 

r 

closed-loop  plant  transfer  function  that  would  raault 
if  k  wara  identically  zero ,  i.a. ,  if  a  constant  control 
law-  kwh*  ware  used.  Under  the  conventional  assumption 
that  tha  plant  relative  degree  is  exactly  known  and,  if 
B„  divides  P,  then  k*  can  be  chosen  [1],  such  that 


. 

A* 


g.,3.. 


ai) 


j*** 


If  the  Relative  Degree  Assumption  is  violated, 

Vm 

can  only  get  as  close  to  — —  as  the  feedback  struc- 

n 

tura  of  tha  controller  allows.  The  first  term  on  the 
right-hand  side  of  eqn.  (9)  results  froa  such  a  consid¬ 
eration.  Note  that  if  eqn.  (11)  were  satisfied,  eqn. 

(9)  reduces  to  the  familiar  error  equation  form  that 
has  appeared  in  the  literature  [8]  for  exact  modeling. 
For  sure  details  tha  reader  is  referred  to  the  litera¬ 
ture  cited  in  this  section  as  mil  as  to  [IS] . 

Figure  2  represents  in  block  diagram  form  the  combina¬ 
tion  of  parameter  adjustment  law  and  error  equations 
described  by  (7)  and  (9) . 

3.  THE  INFINITE  GAIN  OPERATORS 

3.1  Quantitative  Proof  of  Infinite  Gain  for  Operators 
of  CA1 

The  error  system  in  Fig.  2  consists  of  a  forward  linear 
time-invariant  operator  representing  the  nominal  control¬ 
led  plant  complete  with  unmodeled  dynamics,  q*B»  , 

k*A* 

r 

and  a  time-varying  feedback  operator.  It  is  this  feed¬ 
back  operator  which  is  of  immediate  interest.  The 
operator,  reproduced  in  Fig.  3  for  the  case  where  w  is 
a  scalar  and  Tel.  is  parameterized  by  the  function  w(t) 
and  can  be  represented  mathematically  as: 


t 

u(t)' »  Gw(t)  [e(t)  ]■  uQ  +.w(t>  J w(x)e(T)dT 


(12) 


In  order  to  make  the  notion  of  the  gain  of  the  operator 
G  (t) C - j  precise,  we  introduce  the  following  operator 
theoretic  concepts. 

Definition  1:  A  function  fit)  from  (0.™)  to  R  is  said 
to  be  in  if  the  the  truncated  norm 


|f(t> 


Ql 
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(13) 


is  finite  for  all  finite  T. 


Definition  2:  The  gain  of  an  operator  G[f(t)],  which 

maps  functions  in  L,  into  functions  in  L„  is  defined 

as  2*  28 


||G(f<t)l| 


|g| |  «  sup 

f(t,«L2# 

T8(0.») 


f(t) 


(14) 


If  there  is  no  finite  number  satisfying  eqn.  (14) ,  then 
G  is  said  to  have  infinite  gain. 


Theorem  1:  If  w(t)  is  given  by 

w(t)  *  b+c  sins  t 
o 


(IS) 


for  any  positive  constants  b,c,u  ,  tha  operator  of 
eqn.  (12)  has  infinite  gain.  ° 

Proof:  The  proof  consists  of  constructing  a  signal 
e(t) ,  such  that  _ 


G  (elt)  ]  | 


lim 


•  (t> 


(16) 


is  unbounded. 


Let  e(t)  ■  a  sinu  t,  with  a  an  arbitrary  positive 
o 

constant  and  u  the  same  constant  as  in  eqn.  (IS) . 
o 

These  signals  produce: 

w(t)e(t)  *  ab  sinu  t  +  4  ac  -  —  ac  cos2u  t  (17) 

O  <l  i  o 
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Next,  using  standard  norm  inequalities,  we  obtain 
from  eqn.  (19) 


1  l»(t>  |  |^>  |  abett  i  ac2t  sinuQt|  J*  -|  |uj  |*^ 

-1 15s"!  I L  -Ilr£»inw0t||*  '1^-  -  5T-)cosuot||*- 


o  “2 
1 3  abc 


„  »in2uot||  -I,  5S-cos3u0t| 

O  2  O 

>11  J  +  \  »=2t  sinuQt| 


/2 


(20) 


(21) 


with 


«uo|  auj 


-  (22) 


Now 


1 4  abct+4  *=2t  sinu  1 1  |  J 

*  *  O  le- 


?  t a2b2c2  .2  a2c4  .2,2  ..  a2bc3  .2  ,  .  \ 

[~1 -  6  *  T-  *in  “ot+  “ 2~  * 


dt  " 

(23) 
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° 

t  coiaw  t  I 


*sin2u  t 
o 


2  3  ,  _ 

*  *  ■  .  [2u  t  sinu  t-{u  (t)  -2}  cosu  t] 

,  2  o  o  o  o  o 

2“o 


\  / a2b2c2  *  a2c4  ^ 
-  \  12  24  / 


T3  -  K2T2  -  ^T2  -  Ko 


where 


(24) 

<2S) 

(26) 

(27) 

(28) 


Combining  inequalities  (21)  and  (2S)  we  arrive  at: 
.  2 


,*  -Kj*-K  (29) 


Also. 


sin2u  t  dt  <  a  T 
o  — 


(30) 


Therefore, 


'““"s  >[(^  -#)’’-vV. 

|e(t)  I  1 5  _L  a2T 
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and,  therefore,  G^  for  w  aa  in  eqn.  (IS)  haa  infinite 
gain. 


In  addition  to  the  fact  that  the  operator  G^ft)  from 
e(t)  to  u(t)  has  infinite  gain,  the  operator  from 
e(t)  to  fc(t)  in  Fig.  3  also  haa  infinite  gain.  This 
operator  is  described  by: 


infinite  gains  can  arise  froa  any  component  of  the  vec¬ 
tor  w(t) . 

Remark  2:  The  corresponding  operators  Gw  and  de¬ 
fined  for  various  other  adaptive  algorithms  such  as 
the  Marendra,  Lin,  Valavani  [31  and  Horae  [4]  of  the 
model  reference  type,  aa  well  as  the  algorithms  devel¬ 
oped  by  Egardt  [9] ,  which  include  the  self-tuning  re¬ 
gulators,  can  also  be  proven  to  be  infinite  gain 
operators:  sea  Rohrs  [151. 

Remark  3:  Infinite  gain  operators  are  genetically  pres¬ 
ent  in  adaptive  control  and  are  typically  represented 
aa  in  Fig.  4,  where  F(s)  is  a  stable  diagonal  transfer 
function  matrix  and  M  is  (usually)  a  memory less  map. 

D  and  C  are  vectors  of  various  input  and  output  com¬ 
binations,  including  filtered  versions  of  said  signals. 
The  operator  in  Fig.  4  can  also  be  proven  to  be  infinite 
gain  (see  Rohrs  [IS] ) . 

3.2  Two  Mechanisms  of  Instability 

In  this  section,  we  use  the  algorithm  CA1  to  in¬ 
troduce  and  delineate  two  mechanisms  which  may  cause 
unstable  behavior  in  the  adaptive  system  CM,  when  it  is 
implemented  in  the  presence  of  unmodeled  dynamics  and 
excited  by  sinusoidal  reference  inputs  or  by  distur¬ 
bances.  The  arguments  made  for  CA1  are  also  valid  for 
other  classes  of  algorithms  mentioned  in  Remarks  2  and 
3,  mutatis  mutandis.  Since  the  arguments  explaining 
instability  are  somewhat  heuristic  in  nature,  they  are 
verified  by  simulation.  Representative  simulation  re¬ 
sults  are  given  in  Section  4. 

3.2.1  The  Causes  of  possible  Instability 

In  order  to  demonstrate  the  infinite  gain  nature 
of  the  feedback  operator  of  the  error  system  of  CA1  in 
Section  2,  it  is  assumed  that  a  component  of  w(t)  has 
the  fora 

w,  (t)  •  b  +  c  sinu  t  (32) 

i  o 

and  that  the  error  has  the  form 

e(t)  -  a  sinu  t  (33) 

o 

The  arguments  of  Section  2  indicate  that,  if  e(t>  and 
a  component  of  w(t)  have  distinct  sinusoids  at  a  com¬ 
mon  frequency,  the  operator  Gw(t)  sqn*  (12)  and  the 

operator  Hw(fc)  of  eqn.  (31)  will  have  infinite  gains. 

Two  possibilities  for  a(t)  and  w(t)  to  have  the  forms 
of  eqn.  (32)  and  eqn.  (33)  are  now  considered. 

Case  (1) i  If  the  reference  input  consists  of  a  sinu¬ 
soid  and  a  constant,  e.g. 


Hw(t,[.(t)l 


E  + 


o 


w(T)e(T)dT 


(31) 


Theorem  2:  The  operator  with  w(t.)  given  in  eqn. 

(13)  has  infinite  gain. 


r(t)  •  r,  ♦  r,  sinu  t  (34) 

12  O 

where  r,,  and  r  are  constants,  then  the  plant  output 
y(t)  will  contain  a  constant  term  and  a  sinusoid  at 
frequency  u  .  Consequently,  through  eqns.  (2), (3)  and 
(3B),  all  components  of  the  vector  w(t)  will  contain  a 
constant  and  a  sinusoid  of  frequency  u^. 


Proof:  Choose  e(t)  »  a  sinuQt  as  before. 

Then  E(t)  •  Hw(t)[e(t)l  is  given  by  eqn.  (IS). 

Proof  of  Infinite  gain  for  this  operator  then  follows 
in  exactly  analogous  steps  as  in  Theorem  1  and  is, 
therefore,  omitted. 

Remark  1:  loth  operators  G  and  H  will  also  have 
infinite  gain  for  vectors  wYt) ,  since  the  operator 


If  the  controlled  plant  matches  the  model  at  d.c.  but 
not  at  the  frequency  uq,  the  output  error 

eft)  *  y(t)  -  yM<t)  <U> 

will  contain  a  sinusoid  at  frequency  tug.  Thus,  the 

conditions  for  infinite  gain  in  the  feedback  path  of 
Figure  1  have  been  attained. 


Case  (2)  i  If  a  sinusoidal  disturbance,  d(t),  at  frs- 

quancy  a  enter u  tha  plant  output  as  shown  in  Fig.  1, 
the  sinusoid  will  appear  in  w(t)  through  tha  following 
aquation  which  replaces  eqn.  (3)  in  tha  presence  of  an 
output  disturbance 

s1*1 

W  <t>  -  £y(t>*d(tt)];  i-l,2,...,n  (36) 

The  following  equation  replaces  eqn.  (6)  when  an  out¬ 
put  disturbance  is  present 

e(t)  -  y (t)  +  d(t)  -  y  (t)  (37) 

M 

Any  sinusoid  present  in  d(t)  will  also  enter  e(t) 
through  eqn.  (37).  Thus  the  signal  eft)  and  w(t)  will 
contain  sinusoids  of  the  seme  frequency  and  the  op¬ 
erators  Hw^tj  and  Gw<t)  will  display  an  infinite  gain. 

3*2*2  Instability  Due  to  the  Gain  of  the  Operator  a 
of  equation  (12)  — 

The  operator  G  of  eqn.  (12)  is  not  only  an  infinite 
gain  operator  —  but  its  gain  influences  the  system  in 
such  a  manner  as  to  allow  arguments  using  linear  sys¬ 
tem  concepts,  as  outlined  below. 

Assume ,  initially,  that  the  error  signal  is  of  tha  form 
of  eqn.  (33),  i.e.,  a  sinusoid  at  frequency  u  .  Assuan 
also  that  a  component  of  w(t)  is  of  tha  form  Sf  eqn. 
(32),  i.e.,  a  constant  plus  a  sinusoid  at  the  same  fre¬ 
quency  u  as  the  input.  The  output  of  the  infinite 
gain  operator,  of  eqn.  (12),  as  given  by  eqn. 

<19^  consists  of  a  sinusoid  at  frequency  uiq  with  a  gain 

which  increases  linearly  with  time  plus  other  terms  at 
0  radians/sac  (i.e.  d.c.)  and  other  harmonica  of  u  > 

12  0 

1. e.  u(t)  «-ac  t  sinui  t  +  other  terms. 

2  Q 

The  infinite  gain  operator  manifests  its  large  gain 
by  .  roducing  at  tha  output  a  sinusoid  at  the  earn  fre- 
quency,  uq,  aa  tha  input  sinuaoid  but  with  an  amplitude 

increasing  with  time.  By  concentrating  on 

tha  signal  at  fraquancy  and  viewing  tha  operator 

Gw(t)  fig  *  simple  time- Increasing  gain  with  no  phaaa 
shift  at  the  frequency  end  very  smell  gain  at  other 

frequencies .  we  will  be  able  to  come  up  with  a  mech¬ 
anism  for  instability  of  tha  error  system  of  Figure  2, 
where  Gw(e)  consists  of  the  feedback  part  of  the  loop. 

If  the  toward  path,  of  the  error  loop  of  Figure 

*r 

2,  has  less  than  *180*  phaaa  shift  at  tha  frequency  wQ, 

and  if  tha  gain  of  G  ...  were  indeed  small  at  all  other 
w(t) 

frequencies,  then  the  high  gain  of  G  at  u  would 

wlCJ  O 

not  affact  tha  stability  of  tha  error  loop. 

If,  however,  the  torseud  loop,  ^*-**  ,  does  have  180* 

r 

phase  shift  st  u  ,  the  combination  of  this  phase  shift 

with  tha  sign  reversal  will  produce  a  positive  feedback 
loop  around  the  operator  Gw(tj ,  thereby  reinforcing 

the  sinuaoid  at  the  input  of  G„(e)*  The  sinusoid  will 
then  increase  in  amplitude  linearly  with  time,  aa  tha 
gain  of  Gw(t)  grows,  until  ths  combined  gain  of  G^t) 

and  ^7!  exceeds  unity  at  tha  frequency  u  .  At  this 

RJA  O 

point,  the  loop  itself  will  become  unstable  and  all 


signals  will  grow  without  bound  vary  quickly  (aa  tha 
offsets  of  the  unstable  loop  and  continually  growing 
gain  of  Gw(t)  compound.) 

Since  the  infinite  gain  of  G  ...  can  be  achieved  et 
g*B*  *(t) 

any  frequency  u>o»  if  ‘t-,*  has  +180*  phaaa  shift  at 

any  frequency,  the  adaptive  syatam  is  susceptible  to 
instability  froa>  either  a  raferanca  input  or  a  dis¬ 
turbance. 


Thus  tha  importance  of  the  Relative  Degree  Assumption, 


which  assantially  allows  ona  to  aasuma  that  1* 

r 

strictly  positiva  real  is  sean.  Tha  stability  proof 
of  CA1  hinges  on  tha  assumption  that  g«B*  is 

k*A* 

strictly  positive  real  and  that  G^^  Is  passivs,  i.a. 


/ 


G  j  (a(t))a(t)dt  >  0 


(38) 


Both  proparties  of  positive  realness  and  passivity 
are  properties  which  are  independent  of  the  gain  of 
the  operator  involved.  However,  it  la  always  the  case 
that,  due  to  tha  Inavitabla  unmodeled  dynamics,  only  a 
bound  is  known  on  tha  gain  of  ths  plant  at  high  fre¬ 
quencies.  Therefore,  tor  a  large  class  of  unmedelsd 
dynamics  in  tha  plant,  including  all  unaodeled  dynamics 
with  relative  degree  two  or  greater,  the  operator, 

q*B* 

,  will  have  *  180*  phase  shift  et  some  frequency 

and  be  susceptible  to  unstable  behavior  if  subjected 
to  sinusoidal  refaranca  inputs  and/or  disturbances  in 
that  fraquancy  range. 


3.2.3  Instability  Due  to  tha  Gain  of  tha  operator 
Hw  of  Equation  (31) 

In  tha  previous  subsection,  the  situation  was  examined 
where  the  amplitude  of  the  sinusoidal  arror  e(t>  grew 
with  time  due  to  e  positive  feedback  mechanism  in  the 
error  loop,  in  this  subsection,  wa  explore  tha  situa¬ 
tion  where  ths  sinusoidal  error,  s(t),  is  not  at  a 
frequency  where  it  will  grow  due  to  tha  arror  system 
but,  rather,  whan  there  exist  persistent  steady-state 
errors.  Such  a  persistent  arror  could  arise  from 
either  or  both  of  the  two  mechanism  discussed  in 
Section  3.1. 


1)  A  reference  input  with  a  number  of  frequencies  is 
applied  and  tha  controlled  plant  with  unmodalad  dy¬ 
namics  cannot  match  tha  model  in  amplitude  end  phase 
for  all  rafarance  input  frequencies  involved.  This 
will  causa  a  persistant  sinuaoid  in  both  the  arror 
e(t),  through  eqn.  (6),  and  tha  signals  w(t) ,  through 
eqns.  (2)  and  (3) ,  and/or 

(2)  An  output  sinusoidal  disturbance,  d(t) ,  enters 
sa  shown  in  Figure  1,  causing  tha  persiscent  sinusoid 
directly  on  e(t),  through  sqn.  (37),  and  w(t)  through 
aqn.  (36) . 


Aasuma,  that  through  ona  of  tha  above  or  any  other 
mechanism  that  a  component  of  w(t)  contains  a  sinusoid 
at  frequency  u  as  in  sqn.  (32)  and  that  e(t)  contains 
a  sinusoid  of  °tha  sane  frequency.  Then  the  operator 
has  infinite  gain  and  ths  norm  of  ths  output 

signal  of  this  operator,  hit),  inersasss  without  bound. 
Tha  signal,  k(t) ,  will  taka  tha  form  of  sqn.  (18) , 
rapeatsd  here: 


4b  4b  4C 

*  ~ ~  co*‘»ot"43  ,in*V 

00  o 


£<t>  “  K  ♦  4 

o  2 


with  th*  reference  input 
■  r(t)  »  . 3+2.0  ainB.Ot 


(44) 


Proa  tha  aacond  ter*  ana  can  aaa  that  tha  parameters 
of  tha  controller.  defined  in  aqn.  (10),  i.a. , 
k(t)  -  k*  ♦  R(t) ,  will  incraaaa  without  pound. 

If  tharo  ara  any  unmodeled  dynamic*  at  all,  increasing 
tha  aisa  of  tha  nominal  feedback  controllar  para- 
natara  without  bound  will  causa  tha  adaptive  system  to 
became  unatabla.  Indead,  ainea  it  ia  the  gains  of  tha 
Mimi  faadback  loop  that  ara  unbounded,  tha  system 
will  becoom  unatabla  for  a  largo  claaa  of  plant a  in¬ 
cluding  all  thoaa  whoae  relative  degree  ia  three  or 
more,  avan  if  no  uaaodelad  dynaatica  are  praaent. 

4.  SIMULATION  RESULTS 


In  thia  a action  tha  arguments  for  inatability  preaen- 
tad  in  the  previous  sections  are  shown  to  be  valid 
via  simulation. 


Tha  simulation*  ware  generated  using  a  nominally  first 
order  plant  with  a  pair  of  complex  unmodeled  poles, 
described  by 


y(t) 


_2 _ 

(*♦  1) 


229 

(a2+30a+229) 


tu(t)) 


(39) 


and  a  reference  modal 

yM(t)  -  (r(tt)J  (40) 

Tha  simulations  were  all  initializad  with 


ky(0)  --.65  i  kr(0)  -  1.14 


(41) 


This  simulation  demonstrates  that  if  tha  sinusoid  in¬ 
put  is  at  a  frequency  for  which  tha  nominal  controlled 
plant  does  not  generate  a  large  phase  shift  (at 

ui  -8.0,  tha  phase  shift  of  aqn.  (42)  is  -133*),  tha 
o 

algorithm  may  stabilize  despite  tha  high  gain  operator. 

Similar  results  wars  obtained  for  tha  algoritisss  des¬ 
cribed  in  [3,4,6,7,91,  but  ara  not  included  hare  due 
to  space  considerations.  Tha  reader  is  referred  to 
[151  for  a  more  comprehensive  sat  of  simulation  results, 
in  which  instability  occurs  via  both  tha  mechanisms 
described  In  sections  3.2.2  and  3.2.3,  for  sinusoidal 
inputs. 

4.2  Simulations  with  Output  Disturbances 

The  results  in  this  subsection  demonstrate  that  tha 
inatability  mechanism  explained  in  Section  3.2.2  does 
indeed  occur  whan  there  is  an  additive  unknown  output 
disturbance  at  the  wrong  frequency,  entering  tha  system 
as  shown  in  Fig.  1.  In  addition,  the  inatability  mecha¬ 
nism  of  section  3.2.3,  which  will  drive  tha  algorithms 
unstable  whan  there  is  a  sinusoidal  disturbance  at  any 
frequency,  is  also  shown  to  take  place.  The  seam 
numerical  example  is  employed  here  as  well. 


Instability  via  the  Phase  Mechanism  of  Section  3.2.2 

In  this  case,  CA1  was  driven  by  a  constant  reference 
input 

r(t>  -  .3  (45) 


which  yield  a  stable  linearisation  of*  Che  associated 
error  equations.  For  the  parameter  values  of  eqn. 
(41)  one  finds  that 


with  a  very  small  output  additive  disturbance 

d(t>  -  5.59  ■*  10-6  sin  16.lt  (46) 


527 _ 
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The  results  ere  shown  in  Fig.  7,  and  instability  occurs 
(42)  as  predicted.  The  only  surprise  nay  be  the  minuteness 

of  the  disturbance  (  10~b)  which  will  cause  instability. 


Tha  reference  input  signal  was  chosen  baaed  upon  tha 
discussion  of  section  3.2.2 

r(t)  -  .3  ♦  1.85  sin  16.lt  ,  (43) 


the  frequency  16.1  rad/sec.  being  the  frequency  at 
which  the  plant  and  the  transfer  function  in  (42), i. a. 

a*a* 

jL  ,  has  180*  phase  lag.  A  small  d.c.  offset  was  pro* 
r 

vided  so  that  the  linearised  system  would  be  asympto¬ 
tically  stable.  Tha  relatively  large  amplitude,  1.85 
of  the  sinusoid  in  eqn.  (43)  was  chosen  so  that  tha 
unstable  behavior  would  occur  over  a  reasonable  simul¬ 
ation  time.  Tha  adaptation  gains  were  set  equal 
to  unity. 

4.1  Sinusoidal  Reference  Inputs 

Figure  5  shows  the  pient  output  and  parameters  kf(t) 
and  ky(t)  for  tha  conditions  described  so  far.  The 

amplitude  of  the  pleat  output  at  the  critical  frequency 

<u  -16.1  rad/ sec)  end  the  parameters  grow  linearly 
P 

with  time  until  the  loop  gain  of  the  error  system 
beoomas  larger  than  unity.  At  this  point  in  time, 
even  though  the  parameter  values  ere  well  within  the 
region  of  stability  for  the  linearised  system,  highly 
unstable  behavior  results. 

Figure  6  shows  the  results  of  e  simulation,  this  tims 


Instability  via  the  Gain  Incrssss  Mechanism 
of  Section  3.2.3 

Figure  8  shows  the  results  of  a  simulation  of  CA1  that 
waa  generated  with 

r-0. 3 

but  the  disturbance  wee  changed  to 

d(t)  •  8.0x10  6  sin5t  (47) 

cr*B* 

At  u  -5,  of  eqn.  (42)  provides  only  -102*  phase 

shift  to  thi  sinusoidal  error  signal  oC  increasing 
amplitude,  which  is  characteristic  of  Instability  via 
tha  machanism  of  section  3.2.2,  la  not  saan  in  Fig.  8. 
What  ia  saan  is  that  tha  system  becomes  unstable  by 
the  mechanism  of  Section  3.2.3.  While  tha  output  appear 
to  settle  down  to  a  steady  state  sinusoidal  error,  tha 
ky  parameter  drifts  sway  until  the  point  where  the  con¬ 
troller  becomes  unstable.  (Only  the  onset  of  unstable 
behavior  is  shown  in  Figure  8  in  order  to  maintain 
scale) .  We  note  also  that  even  when  the  error  appeared 
settled, its  value  represented  a  large  disturbance  am¬ 
plification  rather  then  disturbance  refection? 

The  most  disconcerting  part  of  this  analysis  is  that 
none  of  tha  systems  analyzed  have  bean  able  to  counter 

thia  parameter  drift  for  a  sinusoidal  disturbance  at 
any  frequency  tried: 


Indeed,  Figure  9  shows  the  result*  of  a  simulation  run 
with  rafaranca  input 

r*0.0  (43) 

and  constant  disturbanca 

d  •  3.0  x  10"6  (49) 

Tha  simulation  rasulta  show  that  ths  output  gain  sattla 
for  a  long  time,  again  with  disturbanca  asgilif ication, 
but  tha  parameter  X  incraasas  in  magnituda  until  ins¬ 
tability  ensues.*  T^us  tha  adaptiva  algorithm  shows  no 
ability  to  act  a van  as  a  regulator  whan  thara  ara  output 
disturbancas. 

5.  COWCLUSIOMS 

In  this  pa par  it  was  shown,  by  analytical  methods  and 
verified  by  simulation  results,  that  existing  adaptive 
algorithms  as  described  in  [1-4,6,7,19],  have  imbedded 
in  thair  adaptation  mechanisms  infinite  gain  operators 
which,  in  tha  presence  of  unmodelad  dynamics,  will 
cause: 

*  instability,  if  tha  reference  input  is  a  high 
frequency  sinusoid 

'disturbance  amplification  and  instability 
if  there  is  a  sinusoidal  output  disturbance 
at  any  frequency  including  d.c. 

'instability,  at  any  frequency  of  reference  inputs 
for  which  there  is  a  non-zero  steady  state  error. 
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While  the  first  problem  can  be  alleviated  by  proper 
limitations  on  tha  class  of  permissible  reference  in¬ 
puts,  the  designer  has  no  control  over  the  additive 
output  disturbances  which  impact  his  system,  or  of  non¬ 
zero  steady-state  errors  that  ara  a  consequence  of 
imperfect  model  matching.  Sinusoidal  disturbancas  and 
inexact  matching  conditions  are  extremely  coesson  in 
practice  and  can  produce  disastrous  instabilities  in 
the  adaptiva  algorithm*  considered. 

Suggested  remedies  in  the  literature  such  as  low  pass 
filtering  of  plant  output  or  error  signal  [26,7,21] 
will  not  work  either. 

It  is  shown  in  [IS]  that  adding  the  filter  to  tha  out¬ 
put  of  the  plant  doe*  nothing  to  change  tha  basic  sta¬ 
bility  problem  as  discussed  in  section  3.2.  It  is 
also  shown  in  [IS]  that  filtering  of  tha  output  error 
merely  results  in  the  destabilizing  input,  being  at  a 
lower  frequency. 

Exactly  analogous  results  were  also  obtained  for  dis¬ 
crete-time  algorithms  as  described  in  [5,17,18,20]  and 
have  bean  reported  in  [15]. 

Finally,  unless  something  is  done  to  eliminate  tha  ad¬ 
verse  reaction  to  disturbancas-at  any  frequancy-and  to 
nonzero  steady-state  errors  in  tha  presence  of  unmodelad 
dynamic*,  tha  existing  adaptive  algorithms  cannot  be 
considered  as  serious  practical  alternative*  to  other 
methods  of  control. 
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Fiqur*  7:  Simulation  of  CA1  with  unmodeled  dynamics, 
r(t)»0.3,  and  d<t)«5.S9xlCT6ainl6.1t. 
(Systam  eventually  becomes  unstable) . 


Figure  9:  Simulation  of  CA1  with  unmodeled  dynamics, 
r(t)»0.0,  and  d(t)-3.0xl0'6. 

(Systam  eventually  becomes  unstable.) 


Figure  8:  Simulation  of  CXI  with  unaodalad  dynamics, 
r (t)“0.3,  and  d(t)«8.0xl0-*«in5.0t. 

(Syatasi  eventually  becomes  unstable) . 


